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We calculate the first passage time distribution for diffusion through a cylindrical pore with 
sticky walls. A particle diffusively explores the interior of the pore through a series of binding 
and unbinding events with the cylinder wall. Through a diagrammatic expansion we obtain first 
passage time statistics for the particle's exit from the pore. Connections between the model and 
nucleocytoplasmic transport in cells are discussed. 
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1. INTRODUCTION 

A current topic of considerable interest in cell biology is determining how cells regulate transport of material 
between the cytoplasm and nucleus [E S H, S Hi- This question is intimately related the study of diffusion in 
confined geometries through the nuclear pore complex (NPC), a roughly cylindrical channel which plays the role of 
gatekeeper to the nucleus. Thousands of NPC's are located on the nuclear envelope, and all material which enters or 
exits the nucleus must pass through individual NPC's. Many of the proteins which make up the NPC are known to 
contain phenylalanine-glycine (FG) repeats, which play an important role in the transport process I, Si- Namely, 
cargo destined for the nucleus form a complex with an importin protein which binds specifically to the FG repeats 
[Qi]. a microscopic model which seeks to describe transport through the NPC should incorporate both the effects of 
geometrical confinement and interactions between the importin-cargo complex and FG repeats lining the NPC. This 
is in contrast to previous studies, which consider the problem within the framework of one dimensional diffusion in 
an effective potential [13, [HI ■ 

In this paper we calculate the first passage time statistics for a particle diffusing in a cylindrical pore with sticky 
walls. The model is motivated by the nucleocytoplasmic transport problem, however we will make a number of 
simplifying assumptions. The first passage problem is formulated in a generic manner which incorporates the effects 
of diffusion in a confined geometry with an attractive interaction at the system boundary. We restrict our attention 
to a cylindrical geometry. Consider a cylinder of radius a and length L, with the axis of the cylinder along the z 
direction (see Fig. [1]). A particle enters the cylinder at z = 0. What is the mean residence time for the particle inside 
the cylinder? What is the distribution of exit times through the cylinder cap located at z = L? 

There are two dimensionless parameters in the problem, the aspect ratio i — —, and y — which is a ratio of 
the characteristic time the particle stays bound to the wall <{, to the time required for the particle to diffuse from the 

2 

center of the cylinder to the wall Here D is the diffusion coefficient of the particle. For a cylinder of length L, 

one could write the mean first passage time through the cyHnder in the scaling form (tc) = jjyg{£,y). In the fimit 
of small £ we have g{i — > 0,y) = 1, since the particle will always reach the end cap before encountering the wall 
and the binding is not important. In the opposite limit of large £ the diffusion is effectively one dimensional. If the 
binding time becomes negligible, we have g{£ — > oo,y — > 0) = 1. We will not pursue a more detailed scaling analysis 
here. Instead, we will demonstrate how the first passage problem can be formulated to treat the case of an arbitrary 
distribution function for the binding time tf, at finite aspect ratio £. To do so, consider drawing all the possible paths 
of the particle through the cylinder. We organize these paths by the number of encounters with the cylinder wall. 
The simplest paths are those for which the particle diffuses to the end of the cylinder without encountering the wall. 
The next set of paths are for particles which diffuse, stick to the wall for some time, detach and diffuse to the end 
of the cylinder. Due to the absorbing boundary conditions on the cylinder wall, the detachment process is modeled 
by radially re-injecting the particle into the cylinder, at which point it resumes diffusing (sec Fig. 3]). As a result. 
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FIG. 1; (Color online). The geometry and boundary conditions for the first passage problem. The cylinder has radius p = a 
and length L along the axial direction z. A particle enters the cylinder at z = and is absorbed at the end cap located at 
z = L. Particles that are absorbed at the side wall will be re-injected into the cylinder (see Fig. |4]), where they resume diffusing 
until encountering the end cap. 

the mean first passage time will depend not only on the distribution of residence times at the cylinder wall, but also 
on the distribution for the particle's radial re-injection. A detailed discussion of this matter is presented in section 
4. We restrict our attention to the case where the particle cannot diffuse along the surface of the wall at p = a. 
This restriction neglects a class of intermittent problems jT2| , where diffusion in the bulk is coupled to diffusion along 
the surface of the cylinder. More complicated paths can be constructed in the same fashion, with more than one 
encounter with the cylinder wall. 

The structure for the paper is the following. In section 2 we solve the diffusion equation inside the cylinder (see 
Fig. [1]), with reflecting boundary conditions at z = 0, and absorbing boundary conditions at z = L and p = a. The 
reflecting boundary condition at z = means that the particle current vanishes on this surface, and all particles 
eventually transit the pore. The choice of absorbing boundary conditions on both the top and bottom surfaces is 
treated in an appendix. The main result is a Fourier-Bessel expansion for the Laplace transform of the Green's function 
(see Eq. (O), from which the particle current and first passage distribution can be determined p^ . In section 3 we 
calculate the first passage distributions to the wall and end cap. As a result we obtain the conditional mean first 
passage time (see Eq. (PH)) ) for particles which traverse the pore without encountering the wall a,t p — a. In section 
4 the formalism is extended to include interactions between the particle and the wall. We solve a recursion relation 
which gives the corrections to the first passage distribution from particle trajectories which include multiple binding 
and unbinding events with the cylinder wall. In section 5 we perform an inverse Laplace transform to determine 
the first passage distribution in the time domain. We discuss how moments of the first passage distribution can be 
obtained for the higher order corrections. 



2. DIFFUSION EQUATION 



Consider the defining relation for the Green's function of the diffusion equation. 




(1) 



Here t = t — t' is the difference between the observation time t and the source time t' , and D is the diffusion coefficient 
for the particle. Rescaling all lengths by L, all times by L^/D, and performing a Laplace transform with respect to 
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the time variable r we have 

(s- v')G(r,5|r') = <5(r-r'). (2) 
Here the Laplace transform of the Green's function is defined as 

G(r,s|r')=/ dTexp(-sr)G(r,r|r'). (3) 
Jo 

We choose to work in cyHndrical coordinates r = {p, z,(f)), so the Laplacian operator takes the following form 



■^^G--— I p— ]+—— + — (4) 
p dp \ dp j p^ d(\P- dz^ 

In the rescaled coordinates, the ends of the cylinder are located at z = and z = 1, and the cylinder wall is located 
at p = l/£ where the dimensionless length scale I = L/a. Eq. ^ can be solved by a Fourier-Bessel expansion (l^. [Tsj 

^ oo oo 

G(r,s|r') = — exp(im((/.-0'))I]-/m(a^mn^P)^(^,^')- (5) 

m— — cxD n— 1 

Here Jm{x) is the Bessel fonction of the first kind of order m, and Xmn is the n*'* root of the equation Jm{x) — 0. 
Inserting Eq. ^ into Eq. ^ yields the following equation for the unknown fimction Z{z, z'), 



<PZ 

dz 



2 ^IrnZ ^ ~A^n{p')5{z ~ z'). (6) 



Above the eigenvalue w,„„ = s + [xmnf)"^ , and the coefficient Amnip') is determined from the Bessel expansion of 
the radial delta function, 

, 21^ f^/^ 5{p - p') 2£^ , 

The solution with the proper symmetry which satisfies the reflecting boundary condition ^ = at z = 0, the 
absorbing boundary condition Z = at z — 1, and has the proper discontinuous first derivative is 

7f /^ Arnn{p')sech{uJ^n) ^ • ^ n ^^ fo\ 

Z[z,z)^ cosh(Wm„z<) smh(Wm„(l - z>)). (8) 

^rnn 



G(r, s|r') = ^ X! exp(im((/) - (j)')) ^ rJ„i{x„uiip)Jm{xmn^p') x (9) 



Here z<, z> is the smaller, larger of z and z'. We arrive at the following expansion for the Green's function 

> prxridmi rh — fh'')'\ > — = 

m— — OO n— i - 

sech(w„„) ^ • r,/ /I ^^ 
cosh(ti;m„z<) smh(w„,.„(l - z>)j. 

With this result, we are in a position to determine the first passage distributions to the wall and end cap. 

3. FIRST PASSAGE DISTRIBUTIONS 

From the Green's function a number of quantities of interest can be calculated. In particular, we require the first 
passage distributions to the wall of the cylinder at p = l/l and to the end cap at z = 1. The current jc{s\r') through 
the end cap can be calculated from the Green's function as 



Us\r') = - d0 / dp 
Jo Jo 



dG(r,s\r') 



dz 



(10) 



2 Jo(a^o«^/oO cosh(won^') (H) 

~^ XOnJiyXOn) 
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FIG. 2: (Color online). The splitting probabilities ttc and tt^, for a particle starting at the origin as a function of the dimensionless 
aspect ratio i. 



Since the Green's function is normalized, the current jclslr') is the first passage distribution to the end cap for a 
particle which starts at r' . The current through the wall at p = l/£ at height z and hence the first passage distribution 
can be calculated as 



j^,(z,s|r') = - 



2-K 



dG(v,s\v') 



dp 



2^2 xonSech(a;orO j^^^^^^^/-) cosh(tJonZ<) sinh(cjo,i(l - Zy)). 
^ Ji(a;o„)wo„ 

71—1 



(12) 
(13) 



The current jw{z, s\v') is the first passage distribution for a particle starting at r' to the wall p — 1/t dX height z. 
By integrating over z we obtain the first passage distribution to the entire wall, 



„1 oo 



a;o„sech(won) j , „ 

:-^Jo{xonip ) X 

{sinh(tJon(l - z'))sinh{u;onz') - cosh(tJo„z')[l ~ cosh(won(l - z'))]}. 



(14) 



These results can be used to obtain the splitting probabilities tt.;, where i G {w,c} indexes the wall and end cap 
respectively. The splitting probability ni is the conditional probability that the particle hits the boundary indexed by 
i. Consider the definition for the Laplace transform of the currents. 



/"OO /"OO 

J,(s|r') = / dr ji(T|r')exp(-sr) = / dTj.,(r|r') 
Jo Jo 

1 - s{t) + 



1 — ST 

2! 



2! 



(15) 
(16) 



Hence we see that the splitting probability tTc = jc(0|r') and all the moments of the first passage distribution 
(t") = ^^"^^ al" ^ ls=o can be obtained by differentiating the Laplace transform of the current. For example. 
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FIG. 3: (Color online). The conditional mean first passage time (tc) (see Eq. (|20|) ') in units of for a particle starting at 

the origin as a function of the dimensionless aspect ratio t. The shaded region covers the range plus and minus one standard 
deviation g — \J (r^) — (tc)^. 

consider a particle starting at the origin p' = z' = 0. The splitting probabilities (see Fig. [2]) ttc and ttw to the end 
cap and wall respectively are (for a particle starting at the origin r' = 0): 

- in\n\ oX^ sech(a:o„^) 
TTc = Jc(0|O) = 2 > — (17) 

11—1 

oo ^ 

= i^(0|0) = 2 V — ^ [1 - sech(xo„£)] (18) 

These splitting probabilities are propertly normalized as evident from an identity of the zeros of the Bessel function 

oo ^ 

71, + ^^ = 2^ ^ = 1. (19) 

n— 1 ^ ^ 

Note that the splitting probabilities, which are independent of the angle </), can also be obtained in a more direct 
manner by solving Laplace's equation inside the cylinder [l3| . For the particle starting at the origin which avoids the 
wall, the average time required to reach the end cap is (see Fig. ^ 



1 tanh(a;o„^)sech(xo„^) ,„„^ 

<^-> = ;r7 — ^r-fi7ZA • (^o) 



Using this result the limiting behavior can be obtained for the mean first passage time. For £ <C 1 the mean first 
passage time originally grows quadratically as (tc) ~ , a result we arrived at in the introductory section using 
scaling arguments. For those paths that do not interact with the wall, in the opposite limit of ^ ;3> 1 we have linear 
dependence on the aspect ratio as (t^) ~ Jj2x^' these two limiting expressions (r^) is given in dimensionful 
units to facilitate comparison with Fig. [3] To quantify the width of the distribution the second moment can also be 
calculated, 



2s _ 1 sech(xon^) 

2^^^'^1^0„^l(^0n) 



l + ^^^%^-2sech2(.o„^) 



(21) 
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4. 



THE HITTING NUMBER EXPANSION 




FIG. 4: (Color online). Graphical depiction of the hitting number expansion. The particle diffuses until it encounters the 
absorbing wall. After some residence time at the wall, the particle is radially re-injected into the cylinder, where it resumes 
diffusing until it exits the cylinder. 

In this section we will calculate the corrections to these results due to interactions between the particle and the 
cylinder wall. We consider a diagrammatic expansion for corrections to the first passage distribution to the end cap 
(see Fig. 2]). As discussed in the introduction, the first such correction is the path for a particle which diffuses, 
becomes attached to the wall for some time, then detaches and resumes normal diffusion. At this point we should 
note that the absorbing boundary conditions at the wall present a difficulty for the calculation. Because the Green's 
function vanishes on the wall, the particle cannot resume diffusing there. In the model, after the absorption phase, 
the particle is displaced radially inward by a finite amount where it resumes diffusing. As a result, including the 
interactions with the wall, the first passage time through the cylinder will depend on two distributions. The first is 
the departure time distribution $(t), which gives the distribution of residence times at the cylinder wall. The final 
result will also depend on the distribution of the starting location for the radial re- injection P{p). Both of these 
distributions will be discussed in more detail below. 

If we consider all diffusive paths, we can calculate the overall first passage distribution /(s|r') to the end cap in the 
following fashion. 



The first term in the series has already been calculated, it is simply the first passage distribution to the wall ]c{s\v'). 
This is a conditional probability distribution, since it provides statistics for particles which avoid the wall during their 
trip to the end cap. 



oo 




(22) 



N=0 



h{s\v')^%{s\v'). 



(23) 



The first order correction can be expressed as 




(24) 



The second order correction is 




(25) 
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FIG. 5: (Color online). First order corrections to the mean first passage time (tc) in units of 

the origin as a function of the dimensionless aspect ratio £. The black curve is the result of Eq. (120^ 
path fo, and the colored curves include the correction /i with different departure rates n (see Eq 
distribution of starting positions P{p) is delta distributed at p, — 0.8/£. 



for a particle starting at 

which considers only the 
l). In the figure, the 



Here we have introduced the Laplace transform $(z, s) of the departure time distribution [16| . The quantity ^{z, T)dT 
is the probability that the particle attached to the wall at height z departs between r and r + dr. In many cases this 
departure distribution is well described by the exponential form p^ . p7j 

$(r) = Kexp(-KT). (26) 

This exponential form is used in Fig. [5l but the general results are valid for any choice of $ which is independent of 
z. The result is (see Fig. [S]) 



/i(.|0) = $(.) )^ (' 2(Jo(xo„^pi)) \ g f2exorn j ^^^^ [sech(u;o„) - sech(c.o™)] . (27) 



Here we have defined 

After detaching from the wall, the particle resumes diffusion through the cylinder. Due to the absorbing boundary 
conditions at p = the particle cannot resume diffusion at the wall since the Green's function vanishes there. For 
this reason we have averaged over the radial position where the particle begins diffusing after detaching from the wall. 
If we let P{p) denote the distribution function for the starting position, then 

n ( w -fo^^ dppP{p)Jo{xonep) 

Jo xo„£p) = -YJ2 • 29) 

J,fdppP{p) 

This model assumption is necessary for the following reasons. After detaching from the wall, a particle which is 
displaced a small distance will encounter the wall again with probability one in an arbitrarily short amount of time. 
This would result in the particle effectively remaining attached to the wall, and exploring only the surface of the 
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FIG. 6: (Color online). First order corrections to the mean first passage time (tc) in units of for a particle starting at 

the origin as a function of the dimensionless aspect ratio £. The departure time distribution is given by Eq. (|26|l with k — 1. 
The black (blue) solid curves are the result with P{p) given by Eq. pO|l with k — 2 and p* = (p* = 0.8/£) respectively. The 
corresponding dashed curves are the results with P{p) = S{p — p,)/p. The inset shows the distribution of starting positions 
P{p) as a function of p in units of [j] as given by Eq. (f30)) with p, — O.S/£ (blue) and p, = (black). 



cylinder. This is contrary to what is observed experimentally in single molecule tracking studies of transport through 
the NPC Q. These studies indicate that the particle spends the majority of its time diffusively exploring the interior 
of the pore. For this reason we have introduced the distribution of starting locations P{p), which allows us to capture 
an essential feature of the transport process. Namely, the particle re-enters the interior of the pore without being 
immediately re-absorbed, and is able to diffuse in the bulk until it encounters another surface binding domain which 
is spatially separated from its last encounter. 

We now indicate how the distribution P{p) could be calculated. One possible scenario is that the absorption is 
modeled as a local potential well at the cylinder wall. The calculation of P(p) can be cast as the problem of determining 
the exit distribution for the noise activated escape from the potential well, which is a well studied problem in the 
theory of stochastic differential equations. For a detailed treatment of this topic for Kramers' escape problem, see 
[Tsl. IT9I [20I. [2]| and references therein. Briefly, the particle's Langevin equation is converted into a phase plane system. 
The domain of attraction at the bottom of the potential well is bounded by a separatrix, which is the locus of points 
for which the particle trajectory is equally likely to return to the well or escape. The exit distribution is then the 
probability distribution function of the position where escaping trajectories hit the separatrix. 

Alternatively, we consider a scenario for determining the distribution P{p) motivated by the structure of the nuclear 
pore complex (NPC). In the physical system, the particle binds to FG repeat regions on polymers within the NPC. 
In this case P{p) is related to the monomer density for the polymers within the NPC. The monomer distribution 
for a polymer brush in a cylindrical pore has recently been studied using simulations and scaling theory [23 |. For 
our present purposes we find that results for the monomer distribution can be approximated using a WeibuU like 
distribution with shape parameter k, scale parameter A, and normalization constant C, 



i_ ^fc-l 
Pip)^C[^] exp 



1 X fc- 

7-P 



(30) 



In Fig. [S]we study the effect of P{p) on the first order corrections to (tc). We compare the results for a delta 



9 



distribution P{p) = &{p — p*) I p with a distribution of the form Eq. ([30|l whith k = 2 and \= {j — p,) (fe^r) ' which 
has its maximum at p^. The quahtative results are very similar. For a delta distribution, shifting the starting location 
from near the wall (p* = 0.8/£) to the center of the pore (p, = 0) results in a slight increase in (tc). As compared 
to the delta distribution, the distribution given by Eq. (|30p results in a small correction. This correction increases 
or decreases (tc) depending on whether walks are weighted more heavily with starting locations p < p* or p > p*, 
respectively. Note that changing P(p) effects not only (tc), but also the probability /i(0|0) for the first order diagram. 
This is to be contrasted with the departure time distribution $(r), which does not effect the diagram probability. 
Changes to the departure time distribution <&(t) have a more pronounced effect on the first order corrections to the 
mean first passage time (see Fig. [5]) than changes to the starting distribution P{p) (see Fig. 15]). 

For a short, fat cylinder, i.e. £ <C 1, the first few terms in the hitting number expansion suffice to provide accurate 
results, since the probability of encountering the wall multiple times before reaching the end cap is quite small. We 
now discuss an iterative procedure to calculate all terms in the hitting number expansion. In principle the calculation 
can be taken to any desired order. We can obtain any order correction recursively from the relation 

fN[s\r')= ( dzi/Ar„i(s|zi)$(zi,s)i^(zi,s|r'). (3f) 
Jo 

In what follows we assume that the departure time distribution is independent of the spatial coordinate z, so that 
$(2:, s) — $(s). In this case the convolution integral can be performed exactly, and we obtain an analytic expression 
for all orders in the hitting number expansion. These considerations lead to the following expression for the N^^^ 
term in the series for a particle which starts at the origin r' = 0, 

( 2(Jo(xo„„^pjv)} \ (jv) 
The + 1 rank tensor A^^'^ (s) can be defined recursively as 

Hence all terms in the series are determined after the first is calculated, which is 

^nino (■5) = ^«ino [sech(tjo„o) - sech(wo„i )] . (34) 

The terms with ni = no must be treated carefully, since direct evaluation gives an indeterminate form. Using 
PHopital's rule we find 

.(1). ^ _ sech(a;o„)tanh(tJo„) 

Using these results, the second order corrections can be obtained directly from the recursion relation Eq. The 
diagonal term must be treated by the limiting procedure described above which gives 

A^nLis) = '"'tJ^r"^ Mi - cosh(2c.o„)) - sinh(2c.o„)) . (36) 
I6wf5„ 

We now have an expression for the Laplace transform of the first passage distribution. In the next section we 
illustrate the procedure for obtaining the first passage distribution back in the time domain, and discuss how the 
moments of the distribution can be obtained by simple differentiation. 



5. THE INVERSE TRANSFORM AND MOMENTS OF THE DISTRIBUTION 



We now have a result for the first passage distribution in Laplace space. Ideally we would like the result back in 
the time domain, which requires taking an inverse Laplace transform, 
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r+ioo 

/iv(r|0) = -^/ d5e^Vjv(s|0). (37) 

ZTTl 

r is a vertical contour in the complex plane chosen so that all the singularities of /jv(s|0) are to the left of it. We 
illustrate the procedure with the first diagram, 

/oNO) = 2f^^^. (38) 

^ XO„Jl(xo„) 

Tl — 1 

To start we note that sech(z) has an infinite number of simple poles at ;z = ^ + tti/c with /c G Z. To exhibit the 
poles explicitly we use the infinite series expansion, 



1=1 4 ' • 



Hence we first evaluate the integral 



— / dssech(u;o„)e"^ (40) 



1 /•r+ioo 
^E(-ir^'(2m-l) / d. 



("mn + S) 



Here we have defined amn = ^^"^ ^ — |- (a;ori^)^- The integrand has a pole at s = — am„. Since all the poles are 
on the negative real axis, we can take F = and close the contour with a semicircular arc parametrized by s = i?e"^ 
with G , In the limit R ^ oo the contribution from the arc vanishes since e'** ~ gi?tcos(</)) ^ q_ -y^g evaluate 
the residue at s = —amn and obtain 

oo 

Imn=7Tj2 (-l)"+'(2m - l)e-«-^ (41) 

m— 1 

Using this result we arrive at an expression for the first passage distribution in the time domain, 

oo ^ CXj 

/o(rlO) ^2nJ2 E (-l)'"+'(2m - l)e— ^ (42) 

So, for the particle which avoids the wall during its journey, the first passage distribution is a sum of exponentials. 
After some transient time, the dominant behavior is exponential with rate constant an = ^ + i^oi^)^- This is 
consistent with the roughly exponential distribution of residence times determined experimentally in single molecule 
tracking studies Q . 

Explicitly performing the inverse transform becomes increasingly difficult for the higher order corrections. As 
discussed earlier, all the moments of the distribution can be obtained by differentiating with respect to s. All of the 
s dependence in Eq. (|32p is contained in $(s) and Al^ti„_^...„j„o(s). The derivatives of ^l^lijv~i...nino (s) can also 
be obtained recursively from Eq. p3|) . 



— (s) = n \—A^^'^'^ (s) - — A(^-i) (s) 

(9s ''""ino njvnAr-l---n2niV / ^ -^nw njv- 1 ■ ■ ■n2no ' 



(43) 



The starting point is 



Once the form of the departure time distribution $ is specified, we obtain moments of arbitrary order analytically by 
computing the derivatives of Eq. ([5^ . 
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6. CONCLUSIONS 

In this paper we calculated the first passage distribution for a particle diffusing through a cylindrical pore with 
sticky walls. We organized the calculation in a diagram series which groups diffusive trajectories by the number of 
encounters with the cylinder wall. By solving the recursion relation Eq. (|3ip . we calculate the contribution to the first 
passage distribution from each diagram. The main result is given by Eq. ([32]), which gives the contribution to the first 
passage distribution in Laplace space from each diagram. This is valid provided that the departure time distribution 
<i>(r) is independent of the spatial coordinate z. Once the form of the departure time distribution has been specified, 
all the moments of the first passage distribution can be calculated using the recursion relation and differentiating (see 
Eq. ((33)) ). There is a single dimensionless aspect ratio 1 = Lja which encodes information about the geometry of the 
problem. We illustrate the procedure for obtaining the first passage distribution in the time domain by performing 
the inverse Laplace transform on the first diagram. Considering only the first diagram in the series, Eq. (|42|) provides 
the probability distribution function for the conditional mean first passage time to the end cap. 

With respect to nucleocytoplasmic transport problem mentioned in the introduction, several comments are in order. 
We have made several assumptions in the formulation of our model. The first passage distribution is obtained for 
the case of unbiased diffusion Because of the reflecting boundary conditions on the top surface, all particles 
eventually transit the pore. In experiments not all of the particles which enter the NPC on the cytoplasmic side will 
successfully make their way to the nucelar side. As a result, the first passage times predicted by the model should 
only be compared with experimental events corresponding to a successful passage through the NPC. An alternative 
choice of boundary conditions, for example both top and bottom surfaces absorbing, can be easily treated within the 
model by making the appropriate changes to the boundary conditions satisfied by Eq. ([5]). The calculation proceeds 
in the same manner as before, although now there will be three splitting probabilities instead of two. The results of 
this calculation are given in an appendix. In the model, the particles bind only to the cylinder wall, which corresponds 
to the EG repeat regions localized on the NPC wall. We hope that this work will provide a basis for the analysis of 
future experimental measurements of the first passage distribution through the nuclear pore complex. 

We acknowledge S. H. Yoshimura, F. Jiilicher, J. S. Bois, and A. V. Klopper for valuable discussions. 
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7. APPENDIX: ABSORBING BOUNDARY CONDITIONS 

In this appendix we derive several results for a modified set of boundary conditions, namely absorbing boundary 
conditions on all surfaces. In the paper we treated the top surface as a reflector. Here we revisit the first passage 
problem for a cylinder where both the top and bottom surfaces of the cylinder are absorbing, as well as the cylinder 
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wall. The logic is the same as in the main text. The unknown function Z{z, z') (see Eq. ([6])) must now satisfy the 
boundary conditions Z = at z — and z = 1. The solution is 

7/ /^ A»«(p')cSCh(cj„„,) 

Z(z,z ) = smh(a;m„z<) smh(cjm„(l - z>)j. (45) 

^rnn 



As a result, the expansion for the Green's function now reads 

> prKn(im.( rh — fh'')'\ > — = 

27r 



I °° °° 21'^ 
G(r, s|r') = ^ X! exp(im((/) - (j)')) ^ r Jm(a;m„^p) J,„(x„„ V) x (46) 



csch(w„„) ^ • I,/ /I ^^ 
smh(w™„z<)smh(wTO„(l - 2;>)). 

There are now three relevant currents. The current through the top surface jt{s\r') at 2 = is 

jMr') - 2 £ j,(^„jp') sinh(a>o„(l - z')). (47) 

n— 1 



The current through the bottom surface Ji>(5|r') at z = 1 is 



ib(s|r') = 2 V '^"'^y^O") Jo(xo„/pO sinh(o.o«^0- (48) 



The current through the wall jw{z, s\r') at p — 1/^ is 



jw{z,s\r') = 2^^ ^Oncsc ( oti) sinh(^on2:<) sinh(a;oK(l - z>))- (49) 

, Jl(a;0n)'^0n 
n— 1 ^ ' 

Integrating over the wall as before we obtain 

- / I '\ o/)2 Xo„CSCh(won) , , „ 

^1 Ji(xo„)w^„ 

{sinh(tJo„(l - z'))[cosh(wonz') - 1] - sinh(wo„z')[l - cosh(u;o«(l - z'))]}. (50) 
As a result we derive the splitting probabilities to the bottom and top of the cylinder 

^fc(r') = J,(0|r') = 2 V ^5^i^^Jo(xo„V) sinh(a;o„^z'), (51) 

n— 1 

^i(r') = it(0|r') = 2 V Jo(xo„V) sinh(a;o„^(l - z' j). (52) 

n— 1 

The sphtting probabihty to the wall is 

7r^(r') = i^(0|r') = 2^^ ^^^^T^.JoMp') x 

{sinh(xo„^(l — z'))[cosh(xo„^z') — 1] — sinh(xo„^z')[l — cosh(xo„^(l — z'))]}- (53) 
One can then average over random initial conditions to obtain the averaged splitting probabilities as follows 

1 /" j3 / / /N 4 ^ csch(xo„£) 



Vol(cylinder) 



/ d3rV,(r') = 1 X] 2!^l^(cosHxoni) - 1). (54) 



For uniform initial conditions as considered above symmetry dictates that iTt = TTt- Averaging over the splitting 
probability to the wall gives 

7f^ = ^ 3 °" (2 - 2cosh(a;o„^) + a;o„^sinh(a;o„^)). (55) 

« „=1 ^On 



The averaged splitting probabilities are properly normalized since 



The mean first passage time to the bottom of the cylinder can be calculated as before. The result is 

-1 dUs\r') 



{n{r')) 



7r(,(r') ds 



s=0 



^ 00 
" TTb(r')£ ^ 



CSCh(xon^) 



Jo{xon(-p') (coth(a;on-^) sinh(xon^2;') — z' cosh(a;on^2')) . 



